home *** CD-ROM | disk | FTP | other *** search
/ Speccy ClassiX 1998 / Speccy ClassiX 98.iso / amiga_system / the_aminet / dev / gcc / libmat.lha / src / mgauss.cc < prev    next >
C/C++ Source or Header  |  1980-01-01  |  749b  |  48 lines
  1. //                MATRIX LIB
  2. //            TOMMY JOHANSSON 1995
  3.  
  4. #include "matrix.h"
  5. Matrix gauss(const Matrix& A)
  6. {
  7.  
  8.  
  9.     int o=1,x,y,i,j;
  10.     float v;
  11.     Matrix C(A);
  12.  
  13.     #ifdef DEBUG 
  14.     printf("Gausseliminerar matris.\n");
  15.     #endif
  16.  
  17.     for(o=1;o<=A.m-1;o++)
  18.     {
  19.         for(x=o+1;x<=A.m;x++)
  20.         {
  21.             for(y=o;y<=A.n;y++)
  22.             {
  23.                 if((y==o)&&(C.koff[o][o]!=0)) 
  24.                     v=C.koff[x][o]/C.koff[o][o];
  25.                 else if(C.koff[o][o]==0)
  26.                 {
  27.                     for(i=o+1;i<=A.m;i++)
  28.                     {
  29.                         if(C.koff[i][o]!=0)
  30.                         {
  31.                             for(j=o;j<=A.n;j++)
  32.                                 C.koff[o][j]+=C.koff[i][j];//lΣgg till rad till rad
  33.                             i=A.m;                           // bryt slingan
  34.                         }
  35.                     }
  36.                 }
  37.                             
  38.                 #ifdef DEBUG
  39.                 printf("v=%f o=%d x=%d y=%d\n",v,o,x,y);
  40.                 #endif
  41.                 
  42.                 C.koff[x][y]-=v*C.koff[o][y];
  43.             }
  44.         }
  45.     }
  46.     return(C);
  47. }
  48.